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It is usually assumed that dark matter is produced during the radiation dominated era. There 
is, however, no direct evidence for radiation domination prior to big-bang nucleosynthesis. Two 
non-standard thermal histories are considered. In one, the low-temperature-reheating scenario, ra- 
diation domination begins as late as ~ 1 MeV, and is preceded by significant entropy generation. 
Thermal axion relic abundances are then suppressed, and cosmological limits to axions are loos- 
ened. For reheating temperatures Trh ^ 35 MeV, the large-scale structure limit to the axion mass 
is lifted. The remaining constraint from the total density of matter is significantly relaxed. Con- 
straints are also relaxed for higher reheating temperatures. In a kination scenario, a more modest 
change to cosmological axion constraints is obtained. Future possible constraints to axions and low- 
temperature reheating from the helium abundance and next-generation large-scale-structure surveys 
are discussed. 

PACS numbers: 14.80.Mz,98.80.-k,95.35.+d 
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I. INTRODUCTION 



The Peccei-Quinn (PQ) solution to the strong-CP 
problem yields the axion, a dark-matter candidate [11, 
0, [l]- If the axion mass TOa ^5 10^^ eV, axions will be 
produced thermally, with cosmological abundance 
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where (?*g^F is the effective number of relativistic degrees 
of freedom when axions freeze out [1, 0, @, [3|- Axions 
with masses in the ~ eV range would contribute to the 
total density in roughly equal proportion to baryons. 

Axions in the ~ eV mass range are relativistic when 
they decouple at Tp = 30 - 50 MeV [1]. Free stream- 
ing then erases density perturbations, suppressing the 
matter power spectrum on scales smaller than the ax- 
ion free-streaming length H, [^, [l3|. Light axions 
would also contribute to the early integrated Sachs- 
Wolfe (ISW) effect [ll|. Data from large-scale struc- 
ture (LSS) surveys and cosmic microwave-background 
(CMB) observations have been used to impose the con- 
straint TOa ^ 1 eV to light hadronic axions [12, [3, [3 • 
These arguments apply to any particle relativistic at 
matter-radiation equality or cosmic microwave back- 
ground (CMB) decoupling, thus imposing the similar 
constraint yi. to^/ i 1 eV to the sum of neutrino masses 

[13, [S mm O, O, HO, HH, m, [H H. 

These constraints rely on abundances computed as- 
suming that radiation domination began earlier than the 
chemical freeze-out of light relics. There is, however, 
no direct evidence for radiation domination prior to big- 
bang nucleosynthesis (BBN) 25]. The transition to ra- 
diation domination may be more gradual than typically 
assumed. In such a modified thermal history, two effects 
may cause relic abundances to change. First, the Hub- 
ble expansion rate scales differently with temperature T 
until radiation domination begins, leading to a different 



freeze-out temperature. Second, entropy may be gener- 
ated, suppressing relic abundances. 

The universe could have reheated to a temperature as 
low as 1 MeV, with standard radiation domination be- 
ginning thereafter [H [13, [H, [H, H HH. This low- 
temperature reheating (LTR) scenario may be modeled 
simply through the entropy-generating decay of a massive 
particle (j) into radiation, with fixed rate F^ and initial 
value Hi of the Hubble parameter. The scalar cf) may 
be the inflaton, oscillating as inflation ends and decaying 
into standard- model particles, or it rnight be a secondary 
scalar, produced during preheating [H, [H, [H, [H, [36| . 
This decay softens the scaling of temperature T with cos- 
mological scale factor a, increasing the Hubble param- 
eter H{T) and leading to earlier freeze-out for certain 
relics. Entropy generation then highly suppresses these 
relic abundances. 

Kination models offer another alternative to the stan- 
dard thermal history, invoking a period of scalar-field 
kinetic-energy dominance [37j . but no entropy produc- 
tion. Without entropy generation, abundances change 
more modestly. 

Past work has shown that cosmological constraints to 
neutrinos, weakly interacting massive particles, and non- 
thermally produced axions are relaxed in LTR [H, [sl, 
[39} . Non-thermally produced axions (TOa ^ 10~^ eV) 
would be produced through coherent oscillations of the 
PQ pseudoscalar d, |3, |43|. In this paper, we obtain 
constraints to thermally-produced hadronic axions in the 
kination and LTR scenarios. While kination modestly 
loosens limits, LTR dramatically changes the cosmologi- 
cally allowed range of axion masses. 

We begin by reviewing these modified thermal histories 
and calculating axion relic abundances. We then gen- 
eralize cosmological constraints to axions, allowing for 
low-temperature reheating and kination. For reheating 
temperatures Tj-h ^ 35 MeV, LSS/CMB Hmits to the ax- 
ion mass are lifted; constraints are also relaxed for higher 



2 



Trh- Constraints from the total matter density are also 
relaxed, but not completely lifted. For r,.h ~ 10 MeV, the 
new constraint is TOa ^1-4 keV, while for Trh — 35 MeV, 
we find that <, 43 eV. If T,h ^ 170 MeV, standard re- 
sults are recovered. After estimating the ability of future 
large-scale-structure surveys to further constrain axion 
masses for a variety of reheating temperatures, we de- 
rive modestly relaxed constraints to axions in the kina- 
tion scenario. We conclude by considering future possible 
constraints to the relativistic energy density of axions in 
a low-temperature reheating model. 



II. TWO NON-STANDARD THERMAL 
HISTORIES: LOW-TEMPERATURE 
REHEATING AND KINATION 

We now review the low-temperature reheating (LTR) 
scenario. We consider the coupled evolution of unstable 
massive particles (j), which drive reheating, and radiation 
R, both in kinetic equilibrium. The relevant distribution 
functions obey a Boltzmann equation with a decay term, 
and may be integrated to yield [H, IH, |4l| : 



1 d [p^a 
^ dT 



1 d[pRa'^ 



r0P0, (2) 



where and pR denote the energy densities in the 
scalar field and radiation, respectively, is the de- 
cay rate of the scalar to radiation, and a is the cos- 
mological scale factor. The evolution of the scale fac- 
tor is given by the Friedmann equation, which is — 
[Stt / {3Mpi)] {p^ + pn) well before matter or vacuum- 
energy domination. The reheating temperature Trh is 
defined by [i,[3llil| 
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where Mpi is the Planck mass and g*,rh is the effective 
number of relativistic degrees of freedom when T = T^h. 
In our calculation of the expansion history in LTR, we 
use calculated using the methods of Refs. la. 1421 ■ as 
tabulated for use in the DarkSUSY package [43|. We 
neglect the axionic contribution to for simplicity and 
assume 3 massless neutrinos. The resulting 10% error 
in 5* leads to a comparable fractional error in the result- 
ing axion relic abundance, and is thus negligible at our 
desired level of accuracy. 

We use dimensionless comoving densities [H, |4l[ : 
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where Tg is the temperature today. At the beginning of 
reheating, (p dominates the energy density and radiation 
is negligible. Thus, as initial conditions, we use pr = 
and p^ = [3/ (Stt)] M^^Hf, where Hi is the initial value of 
the Hubble parameter. The two physical free parameters 



in this model are Tih and Hi. The temperature is related 
to the radiation energy density by [6j] 



T = 
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pR ■ 
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We numerically integrate Eqs. ^ to obtain the depen- 
dence of T on a, and the results are shown in Fig. [TJ As 
the scalar begins to decay, the temperature rises sharply 
to a maximum at = (8/3)^^^ aj, where aj is the initial 



value of the scale factor, and then falls as T 
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This shallow scaling of temperature with scale factor re- 
sults from the continual dumping of scalar-field energy 
into radiation, and yields an unusually steep dependence 
of scale factor on temperature. As shown in Figs. [T] and 
[21 when the comoving radiation energy density R over- 
takes $ near T ^ Tih, the epoch of radiation domination 
begins, with the usual T cx scaling. 

Well before reheating concludes, $ is constant and 
PR ^ P4>- If flm ^ a < a (Trh), an approximate solu- 
tion of Eqs. © for T(a) is then [H, EH 
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where 5*,rii = 9* (21h) and Tmax is the maximum temper- 
ature obtained (see Fig. [1]). Here Trh,McV and Hi ^v are 
the reheating temperature and initial value of the Hubble 
parameter, in units of MeV and eV, respectively. 

During reheating, the Hubble parameter is given by 
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^^'9l (T) 
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At a given temperature, the universe thus expands faster 
during reheating than it would during radiation domina- 
tion, and the equilibrium condition T = n (av) > H is 
harder to establish and maintain. Relics with freeze-out 
temperature Tp > Tmax will thus have highly suppressed 
abundances because they never come into chemical equi- 
librium. Relics with Trh ^ Tp <, T^ax come into chemical 
equilibrium, but then freeze out before reheating com- 
pletes. Their abundances are then reduced by entropy 
production during reheating. In either case, species with 
Tp <; Trh have highly suppressed relic abundances. 

Less radical changes to abundances follow in kina- 
tion scenarios. During epochs dominated by the ki- 
netic energy of a scalar field, the energy density p scales 
according to p cx 0"^ [371. Thus H (T) cx T^ or 
H ~ i/rad (T) (T/Tkin), whcrc Frad (T) is the standard 
radiation-dominated H (T) and T^in denotes the transi- 
tion temperature from kination to radiation domination. 
Kination yields relic freeze-out temperatures somewhere 
between the standard and LTR values. There is, how- 
ever, no entropy generation during kination, leading to 
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FIG. 1: This plot shows the evolution of temperature with 
scale factor in a low-temperature reheating (LTR) cosmology 
with Trh = 20 MeV and 3 different initial values for the Hub- 
ble parameter Hj. After a rapid rise due to </> decay, T oc a^''''* 
until T ~ Trh, after which radiation domination begins, and 
T oc a~^. The small bump near T ~ 200 MeV results from a 
jump in near the QCD phase transition. 
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FIG. 2: The top panel shows the Hubble parameter as a func- 
tion of temperature in an LTR cosmology with Trh ~ 20 MeV 
and Hi ^ 1 MeV. Initially H oc T^p, (T), but at tempera- 
tures cooler than T ^ Trh, H oc ^ (r).The second panel 
shows the comoving radiation energy density R = pnT^'^a"' 
and scalar energy density "1> = P4>T^'^T~^ o?" as a function of 
temperature. At T ~ T^h, R flattens out to a constant and $ 
drops ofi' to zero, indicating the conclusion of reheating. 



a less dramatic change in relic abundances. Note that 
these conclusions are rather general, as we have not re- 
lied on any detailed properties of the kination model, but 
only on the scaling (T) cx [H. 



III. AXION PRODUCTION IN 
NON-STANDARD THERMAL HISTORIES 

Axions with > 10^^ eV are thermally produced 
in the standard radiation-dominated cosmology. We now 
show that in LTR models, these axions have suppressed 
abundances. We consider standard hadronic axions, 
which do not couple to standard-model quarks and lep- 
tons at tree level but do have higher-order couplings to pi- 
ous and photons [H, . We do not consider flaton mod- 
els, or other scenarios in which PQ symmetry breaking 
is related to supersymmetry breaking [H, [13, IH, [H, [sO] ■ 
For temperatures T <, 150 MeV, the dominant channels 
for axion production are 7r+ + tt^ a + tt^ , 7r+ + tt'^ — s- 
TT"*" -|- a, and 7r~ + ir'^ ^ a + tt^. Axion scattering rates 
are suppressed relative to particle-number-changing in- 
teractions by factors of T"^ / fpq and thus decouple very 
early. Thus, axions stay in kinetic equilibrium because 
ofTT + TT ^ a + TT, and kinetically decouple when they 
chemically freeze out. 

Nucleonic channels are negligible at these tempera- 
tures. If pious are in chemical equilibrium and Bose en- 
hancement can be neglected, the axion production rate 
is [4, 6, 12, 5l 
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s = 2 [ml + {yiy2 - XlX2^J)\ , (8) 



where — pi /T is the dimensionless pion momenta, yi = 
\/xf + m^/T^ is the dimensionless pion energy, f (yi) = 
1/ [exp(y,i) — 1] is the pion distribution function, Cajr is 
the dimensionless axion-pion coupling constant, and 
is the Riemann ^-function. The energy scale /a = fpq/N 
is the PQ scale, normalized by the PQ color anomaly N, 
and m^r — 135 MeV is the mass of a neutral pion [s^ . 
The PQ scale can be expressed in terms of the axion mass 
fH: 



/a ■ 



1 + r TTis, 



(9) 



Here r = m-a/m^ ~ 0.56 is the up/down quark mass ratio 
and /tt ~ 93 MeV is the pion decay constant [6, 52]. 

Evaluating Eq. ([5]) for P and numerically solving 
Eq. ([2]) for H (T), we estimate the axion freeze-out tem- 
perature using the condition F (Tp) ~ H{Tf). As the 
reheating temperature is lowered, axions freeze out at 
higher temperatures due to the higher value of H, as 
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shown in Fig. [3l As the reheating temperature is in- 
creased, the T oc a~^/^ epoch becomes increasingly ir- 
relevant, and the freeze-out temperature of the axion 
asymptotes to its standard radiation-dominated value. 
Examining Eq. ([8]), we see that F oc oc m^, so higher- 
mass axions keep up with the Hubble expansion for longer 
and generally decouple at lower temperatures. Thus, for 
higher rria, a more radical change to the thermal history 
(even lower Trh), is needed to drive T-p to a fixed higher 
value. 
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3; Freeze-out temperature of the reactions tt -|- tt" <-> 
TT" -I- a, TT^ -f tt" ^ Tr"*" -I- a, and 7r~ -I- tt" ^ vr" -I- a, shown 
as a function of the reheating temperature Trh, for 4 different 
axion masses. More massive axions are coupled more strongly, 
leading to later freeze-out than for lighter axions. 

As axions are spin-0 relativistic bosons, their num- 
ber density at freeze-out is ria (7f ) = C(3)TF/7r^. If 
we assume that axion production ceases at freeze-out, 
the density of axions at any subsequent time is just 
JT-a (Tf) [ap/oo]^, where a-p is the value of the cosmologi- 
cal scale factor at axion freeze-out. The reheating time 
scale, trh — 1 /r, is much shorter than the Hubble time for 
T < Trh, and so it is a good approximation to treat the 
break between the T oc a~^/^ and T oc epochs as in- 
stantaneous at T = Trh. Doing so, we apply Eq. ^ prior 
to the completion of reheating and a oc T^^g'^^^^ (T) af- 
terwards, to obtain 



where rria is the axion mass in units of eV. 

Low reheating temperatures drive up the freeze-out 
temperature. When Trh ^ 7f, the present mass density 
in axions is severely suppressed, because of the sharper 
dependence of the scale factor a on T during reheating. 
This is a result of entropy generation. Using the numeri- 
cal solution for a (T) from Sec. |TT1 we obtain Oa, account- 
ing for the smooth transition between the T oc a~^^^ and 
T oc a^^ regimes. In Fig. |4l we show f^a normalized by 
its standard value, fl^, as a function of Trh. At reheat 
temperatures just a factor of a few below the usual axion 
freeze-out temperature for a given axion mass, the axion 
abundance is suppressed by a factor of 10"'' — 10~^. For 
Trh S> Tp, the axion abundance asymptotes towards its 
standard value. 

In the case of kination, axion freeze-out temperatures 
are still raised, but there is no additional entropy produc- 
tion. Axion abundances are given by Eq. ([1]), but with 
the higher p values appropriate at higher values of 
Tp. 

For the LTR case, our results do not depend on the 
initial value Hj of the Hubble parameter. As seen in 
Fig. [U changes to Hi determine the moment of the fast 
rise to Tmax, but have little influence on the expansion 
history for T < Tmax- For convenience, we choose Hi = 
1 MeV for our calculations, corresponding to Tmax — 
20 GeV. 

Our calculation is valid only if axions are produced in 
equilibrium by thermal pions. This requirement imposes 
the restriction Trh /5 10 MeV. Outside this range, our 
assumptions break down. For sufficiently low values of 
rria and Tj-h, pionic cross sections lead to Tp ^ 200 MeV, 
earlier than the quark-hadron phase transition. The ab- 
sence of hadrons then necessitates the use of quark-axion 
production cross sections. 

Furthermore, for Trh ^ 10 MeV, pions will decay be- 
fore they can come into equilibrium. In both cases, ax- 
ion abundances are suppressed relative to our calcula- 
tion. For axion masses saturating our upper limits and 
Trh <; 10 MeV, we have checked that we are well within 
the equilibrium regime. We restrict ourselves to this 
range, noting that for Trh ^ 10 MeV, more suppressed 
abundances will lead to an even more dramatic relax- 
ation to cosmological axion limits. Finally, coherent os- 
cillations of the axion field pro duce a condensate that 
behaves as cold dark matter [38| , but the resulting addi- 
tional abundance is negligible for ma ^ 10"^ eV at all 
values of Trh under consideration here [ssl . 



IV. CONSTRAINTS TO AXIONS 
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Most constraints to the axion mass are obtained from 
its two-photon interaction. This interaction is parame- 
terized by a coupling constant QaTT; given by ^, i6|, ,5^ . 
y , |56j, l57|, IM l52|, led, l61|, 16^, 1^ 
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FIG. 4: Axion abundance f2a normalized by its standard value 
fia, shown as a function of the reheating temperature Trh. 
Curves are shown for 4 different axion masses. More massive 
axions freeze out later. As a result, their density is less diluted 
by entropy production during the reheating epoch, and so 
they have higher relic densities. 



where ^ is a model-dependent parameter and a is the 
fine-structure constant. The tightest constraint to gs^-yy, 



gaj-y ^ 0.6 X 10"^" GeV"\ 



(12) 



comes from the helium burning lifetime of stars in star 
clusters [HI, |63| . The parameter ^ is given by ^5^, (55|, |56|, 



E 2 (4 + r) 
N ~ 3 (1+7) 



(13) 



Here E and N are weighted sums of the electric and PQ 
charges of fermions that couple to axions. In existing 
axion models, < E/N < 8/3 ■ |65| ■ while r is poorly 
constrained and lies in the range [65l. [66l. l&A [11] : 



0.2 < r < 0.8. 



(14) 



As a result, for any axion model in the range 1.93 ^ 
E/N < 2.39, there are experimentally allowed r values 
for which ga-yy vanishes (see Eg. I13|) . and so constraints 
to axions from star clusters, helioscope, RF cavity, and 
telescope searches may all be lifted [3, [63 . Such a can- 
cellation is fine tuned, but even for other values of E/N, 
constraints to the axion mass are loosened. 

In contrast, the hadronic couplings do not vanish for 
any experimentally allowed r values Q. Axion searches 
based on these couplings are underway, and have already 
placed new upper limits to the axion mass in the keV 



range [69| (and references therein). These couplings also 
determine the cosmological abundance of axions, and so 
useful constraints may be obtained from cosmology. 

Although the hadronic coupling determines the relic 
abundance of axions, ^ will determine the lifetime of the 
axion, which may have implications for cosmological con- 
straints. For the high axion masses allowed by our new 
constraints and certain values of r and E/N, the decay 
a — > 77 may no longer be negligible on cosmological time 
scales. Such an axion would be completely unconstrained 
by limits to ilah^ from the total matter density or the hot 
dark matter mass fraction. In the following calculation, 
we neglect axion decay. Consistent with a vanishing two- 
photon coupling for E/N = 2, we use the value r = 0.50. 
We have verified that our results for flah^ vary by only 
5% for variations in r of about 20%, as the dependence 
of the axion production rate and Tp on r is weak (see 
Eq. [5]) compared with the dependence on the e"™"/^ 
Boltzmann factor. 



A. Constraints to the axion mass from Q^h 

In the standard cosmology, a conservative constraint 
is obtained by requiring that axions not exceed the total 
matter density of Q^h^ ~ 0.135 [l^l, yielding the limit 
TOa ^ 22 eV, using a concordance value for the dimen- 
sionless Hubble parameter h = 0.73. In LTR scenar- 
ios, axion abundances are highly suppressed, as shown in 
Eq. PH)) and Fig. [H Mass constraints to thermal axions 
from cosmology are thus considerably relaxed. 

To obtain abundances in the LTR scenario, we apply 
the numerical freeze-out and abundance calculation of 
Section IIIII Axion mass limits resulting from the re- 
quirement that axions not exceed the total dark matter 
density are shown by the dot-dashed hashed region in 
Fig. O If Trh ^ 40 MeV, constraints are considerably 
relaxed. For example, if Trh — 10 MeV, the axion mass 
constraint is ma < 1.4 keV. When Tj-h ^ 95 MeV, we 
obtain ma ^ 22 eV, equal to the standard radiation- 
dominated result. As already discussed, abundances are 
only slightly changed in the case of kination, so fia^^ 
imposes the constraint ma ^ 26 cV if Tkin — 10 MeV. 



B. Constraints to the axion mass from CMB/LSS 

data 

Axions will free stream at early times, decreasing the 
matter power spectrum on length scales smaller than 
the comoving free-streaming scale, evaluated at matter- 
radiation equality. For sufRciently low masses, axions will 
also contribute to an enhanced early ISW effect [11] in the 
CMB. This suppression is given by AP/P ~ Sfl^/ilm if 
r^a ^ p^ll3l.F70| (and references therein). The mat- 
ter power spectrum thus imposes a constraint to flg.h'^. 

First we review the constraints imposed by these effects 
in a standard thermal history. In this case, both fl^h^ 
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and the free-streaming scale, Afs, depend only on ma in 
a hadronic axion model. Using Sloan Digital Sky Survey 
(SDSS) measurements of the galaxy power spectrum [7i| 
and Wilkinson Microwave Anisotropy Probe (WMAP) 
[T^ I f^^-year measurements of the CMB angular power 
spectrum, Refs. [H, [13, llSl derived hmits of ma ^ 1 eV. 
Axions in the mass range of interest (ma ^ 1 eV) are 
non-relativistic at photon-baryon decoupling, and so this 
constraint essentially comes from measurements of the 
galaxy power spectrum [l^. As a result, we do not ex- 
pect that an analysis including more recent WMAP re- 
sults would make a substantial difference in the allowed 
axion parameter space. In order to lift this constraint, 
the relationship between ma and fls^h^ or Afs must be 
changed. 

If Tp is allowed to vary freely, the constraints may be 
relaxed. In particular, using Eq. ([T]) we can see that in- 
creasing (7*s,F (and thus Tp) will decrease flah^- Fur- 
thermore, if the free-streaming length of the axion is 
less than the smallest length scale on which the lin- 
ear power-spectrum may be reliably measured (Amin — 
40 Mpc), its effect on the matter power spectrum is 
not observable [l^, [l^ [T^I • The comoving free-streaming 
length scale at matter-radiation equality^ @, 
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FIG. 5: Upper limits to the hadronic axion mass from cosmol- 
ogy, allowing the possibility of a low-temperature-reheating 
scenario. The dot-dashed hatched region shows the region 
excluded by the constraint Qi^h^ < 0.135 as a function of 
reheating temperature Trh- The solid hatched region shows 
the additional part of axion parameter space excluded by 
WMAPl/SDSS data. At low reheating temperatures, upper 
limits to the axion mass are loosened. For Trh <; 170 MeV, 
the usual constraints are recovered. 



is set by the ratio between the axion and neutrino tem- 
peratures. 



^ « (10.75/5,3,p)l/^ 



(16) 



so that if 5*3, F <^ 87 (i.e., if axions freeze out considerably 
before neutrinos), the constraint to axion masses is lifted 

m- 

In the case of a modified thermal history, the rela- 
tionship between Tp and ma acquires dependence on an 
additional parameter (Tj-h, in the case of LTR, or Tkin, in 
the case of kination) thus allowing us to loosen the con- 
straints. At a series of values of 5*3, p, Refs. [H, [H, [t^ 
determine the maximum values of fla.h^ consistent with 
WMAP measurements of CMB power spectra and SDSS 
measurements of the galaxy power spectrum. We be- 
gin by mapping these contours, from Fig. 5b in Ref. [l^ 
(which do not include constraints from the Lyman-a for- 
est), into the (ytah^,X{s) plane. For a fixed ma or fia, 
Afs scales monotonically with (?*s.f, and thus serves as a 
proxy for 5,3, f- 



^ This differs from the expression in Refs. [lol . l7Cll | because we as- 
sume, as is the case in our parameter region of interest, that 
m > Teq, the temperature at matter-radiation equality. In 
Ref. [id . m < Tcq is assumed. 



In the domain 10 MeV < T,.h < 250 MeV and 
0.01 eV ^ ma ^ 10 keV, we calculate fl^ (Tlh, Wa) /i^ for 
hadronic axions in LTR, using the full numerical solution 
described in Sec. IIIII We also calculate Afs (Trh, ma). 
Since axions freeze out while relativistic, their energy 
will redshift as E cx a"^. They will have temperature 
Ta = Tpap/a. Meanwhile, the temperature of the cou- 
pled radiation redshifts as T oc a~^/® until radiation 
domination begins. Thus entropy generation modifies 
the relationship between the axion and neutrino temper- 
atures to 



Ta 



1/3 



5/3 



5*,rh5*s,0 
.9*,F3*s,i'h 



1/3 



(17) 



if Tp > Trh. To obtain all of our constraints we use the 
more precise scaling accounting for the smooth transition 
between the T cx a~^/^ and T cx a^^ regimes. The domi- 
nant change to the free-streaming length comes from the 
modified axion temperature, while the modified expan- 
sion rate itself induces negligible fractional changes of 
order TNR/Trh, where Tnr is the cosmic temperature at 
which the axion goes non-relativistic. 

For each pair (Trh, ma), we calculate ila.h^ and Afs 
to trace out the region forbidden with 95% confidence. 
When fia (Trh, ma) /i^ > 0.014, outside the domain of 
Ref. [T^I, we extrapolate, assuming that the 95% con- 
tour asymptotes to a line of constant axion free-streaming 
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wavelength Afs = 40 Mpc. Such a trend is noted in 
Ref. [701], and at the maximum value of ^Ig^h^ of the con- 
tour obtained from Ref. , the maximum allowed free- 
streaming length is consistent with our assumed asymp- 
tote. 

We obtain the upper limit to the axion mass as a func- 
tion of Tih, shown in Fig. [5] Existing LSS/CMB con- 
straints are severely relaxed in the LTR scenario, and 
lifted completely for T,i, ^ 35 MeV. For T,^ <, 35 MeV, 
Afs < 40 Mpc for Wa ^ 1 eV, and so the axion 
mass is unconstrained. It will still be subject to phase- 
space constraints if it saturates the bound fia/i^ ^ 0.135, 
and is to compose all the dark matter in galactic halos 
[zl, 111, llll. At high reheating temperatures, the con- 
straint from LSS/CMB data {n^h'^ <, 0.006) supercedes 
the constraint fl^h^ ^ flj^h^- 

The narrow allowed region between the LSS/CMB and 
total matter density constraints in Fig. \5\ (45 MeV ^ 
Trh ^ 55 MeV) may be simply understood. Axions in 
this narrow window are cold and massive enough to evade 
large-scale structure constraints (i.e., Afs < Amin), and 
dilute enough to evade constraints from the total matter 
density. We note that the CMB/LSS limits asymptote to 
their standard value of Wa ^ 1.4 eV for Trh 170 MeV. 

Future instruments, such as the Large Synoptic Sur- 
vey Telescope (LSST), will measure the matter power- 
spectrum with unprecedented precision {AP/P ~ 0.01) 
[761 177|. This order of magnitude improvement over past 
work [78l [79| will improve the constraint to flah^ by an 
order of magnitude, resulting in the improved sensitivity 
to axion masses and reheating temperatures shown by 
the dotted line in Fig. [51 To estimate possible constraints 
to axions from LSST measurements of the power spec- 
trum, we recalculated our limits using the approximate 
scaling AP/P ~ -8fJa/f^m, assuming AP/P - 10"^ for 
A > 40 h-^ Mpc. 

We also estimate the possible improvement offered 
by including information on smaller scales (Amin ~ 
1 Mpc), as may be obtained from measurements of 
the Lyman-a flux power spectrum [s^l, also shown in 
Fig. [5] We include this effect by replacing Amin with this 
lower minimum length scale. This is indicated by the 
dashed line in Fig. [Sj We can see that more massive ax- 
ions are probed because of information on smaller length 
scales, as are lower reheating temperatures. 

In the case of kination, a much less severe relaxation 
of limits to axions is obtained. As there is no entropy 
generation in the kination case, the abundance and tem- 
perature of the axion are still given by Eqs. (|T]) and 
with the value of 5*3, f appropriate at the new freeze- 
out temperature. In the range of parameter space ex- 
plored, 10 MeV <, Tf <, 100 MeV, and so the vari- 
ation in g*s.F as a result of kination is ~ 60%. For 
Titin — 10 MeV, the new allowed regions are nia ^ 3.2 eV 
and 17 eV ^ ma ^ 26 eV. These conclusions apply to 
any non-entropy-generating scenario in which H (x 
at some early epoch, and not only to kination [25|. If 
Tih ^110 MeV, standard results are recovered. 
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FIG. 6: Estimated improvement in the accessible axion pa- 
rameter space from including more precise measurements of 
the matter power spectrum (region bounded by the dot- 
ted line), corresponding to LSST ffdl. [77j. or from measure- 
ments of clustering on smaller length scales, corresponding to 
Lyman-a forest measurements (region bounded by the dashed 
line) [131 ■ The hatched region indicates the parameter space 
excluded by WMAPl/SDSS measurements. 



V. AXIONS AS RELATIVISTIC DEGREES OF 
FREEDOM AT EARLY TIMES 

Future limits to axions in the standard radiation- 
dominated and LTR thermal histories may follow from 
constraints to their contribution to the energy density 
in relativistic particles at T ~ 1 MeV. Axions are rela- 
tivistic spin-0 bosons, and so pa — (7r^/30) Tp (ap/o)"* = 

(7r^/30) Tp (ap/ayh)* (orh/a)^ H- We can express the to- 
tal relativistic energy density in terms of an effective neu- 
trino number 



n: 



iff _ 



Pa + Piy 
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4/3 



X 3 X 



15 



(18) 



Treating the transition between the T oc a^^^^ and T oc 
a^^ epoch as instantaneous, we solve for the photon and 
axion temperatures, and then obtain 



*(2/) 



= 3 - 



|(f)'/'*(TF/Trh), 



5*s,F 



-4/3 



- 1 



-4/3 



if y > 1, 

if y < 1. 



(19) 
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FIG. 7: Total effective neutrino number jV^ff-m^^" for axions 
with masses saturating the tightest bound on axion masses 
from Fig. [5] The requisite higher temperatures lead to earlier 
axion freeze out and lower Nf The thick black line 
indicates the anticipated sensitivity of CMBPol fs^] to A*'^'^ 
through the primordial helium abundance. 



translate measurements of Yp to constraints on ma and 
Tih we use the scaling relation [s^ 



ANf = 



y{(6.25Ayp + lf -l} 



(20) 



Constraints to N^^ from direct measurements of Yp, in- 
cluding a determination of from CMB observations, 
lead to the 68% confidence level upper limit of Nf < 
3.85 [ii,!!!!!^. From Fig. [7] and Eq. we see that 

this bound cannot constrain rria or Tih- If future mea- 
surements reduce systematic errors, constraints to Trh 
will be obtained for the lighter-mass axions. 

Constraints to ma and Tih may also follow from in- 
direct CMB measurements of Yp. The presence of ^He 
affects CMB anisotromes by changing the ionization his- 
tory of the universe [83| • The Planck satellite is expected 
to reach Al^p — 0.013, yielding a sensitivity of < 
4.04, while CMBPol (a proposed future CMB polariza- 
tion experiment) is expected to approach A Yr, = 0.0039, 
leading to the sensitivity limit Nf < 3.30 [HM 113, IM| ■ 
As shown in Fig. [71 for T^h ^ 15 MeV, such measure- 
ments of Yp may impose more stringent limits on the 
axion mass. Also, if axions with mass in the eV range 
are directly detected, Yp might impose a surprising upper 

limit to [el,!!!. 



For sufficiently high masses, the axionic contribution sat- 
urates to SNf = 4/7 at high reheating temperatures [I]- 
In Fig. El we show Nf''^'^''{T,h), the effective neutrino 
number evaluated at the axion mass which saturates the 
LSS/CMB bounds, for T^h ^ 35 MeV, or saturates the 
constraint ^lah^ ^ 0.135 for lower Tj-h. The behavior of 
the curve may be readily understood. As can be seen 
from Fig. |5l as we increase Tj-h, the maximum allowed 
ma decreases. For T^h ^ 20 MeV, even though the max- 
imum allowed ma is large (which corresponds to a lower 
Tp, since F cx ma), the amount of entropy production 
between Tp and Trh leads to a small axionic contribution 
to Ncff. As Tih increases, the interval between freeze-out 
and reheating decreases. This lessens the impact of en- 
tropy generation, and leads to the rise in A'cff. Finally, 
for Trh <J 20 MeV, the impact of entropy generation is 
nearly negligible, and Ngg falls as the maximum allowed 
value of ma decreases, as in the standard case (due to 
earlier freeze-out). 

A comparison between the abundance of "^He and the 
predicted abundance from BBN places constraints to the 
radiative content of the Universe at T ~ 1 MeV [8l| ; this 
can be stated as a constraint to Nf . At early times, ax- 
ions will contribute to the total relativistic energy density 
(through Nf ), and thus constraints to ^He abundances 
can be turned into constraints on ma and T^h, as shown 
in Fig. m 

In terms of the baryon- number density Ub, we write the 
primordial "^He abundance as Yp = Annc/n-b- In order to 



VI. CONCLUSIONS 

The lack of direct evidence for radiation domination 
at temperatures hotter than 1 MeV has motivated the 
introduction of kination, low-temperature reheating, and 
other scenarios for an altered pre-BBN expansion his- 
tory. In the case of kination, the change in axion abun- 
dances and thus cosmological constraints is modest. Low- 
temperature reheating will suppress the abundance of 
thermally-produced hadronic axions, once the reheating 
temperature T^h ''^ 50 MeV. This is rather intuitive once 
we recall that the axion freeze-out temperature in a radia- 
tion dominated cosmology is ~ 50 MeV. If the reheating 
temperature crosses this threshold, axion densities are 
severely reduced by dramatic entropy production during 
reheating. 

Total density, large-scale structure, and microwave 
background constraints to axions are all severely loos- 
ened as a result, possibly pushing the the axion mass 
window to very high values; for T^h — 10 MeV, the new 
constraint is ma < 1.4 keV. For T^h 170 MeV, stan- 
dard radiation dominated results are recovered. The in- 
clusion of information on smaller scales will probe higher 
axion masses and lower reheat temperatures. More pre- 
cise measurements of the matter power spectrum on all 
scales will probe lower axion masses. Kination also re- 
laxes constraints to axions, though much less markedly. 
Future probes of primordial helium abundance will ei- 
ther lead to further constraints on axion properties, or, 
if axions are directly detected, provide a new view into 
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the thermal history of the universe during the epoch 
10 MeV <,T <, 170 MeV. 
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